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摘要 光 钟 的 频率 稳定 度 和 不 确定 度 达到 了 10-18 量 级 , 使 其 有 望 成 为 下 一 代 的 时 间 频 率 标准 , 并 可 能 用 来 重 
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Vondrak-Cepek 组 合 滤波 算法 应 用 在 光 钟 与 氨 钟 联 
合计 算 的 时 间 尺 度 上 , 尝试 解决 不 连续 的 光 钟 数据 
参与 时 间 尺 度 计 算 的 难点 问题 . 

本 文 首先 对 氧 钟 和 光 钟 数据 进行 了 仿真 , 采用 
ALGOS 算 法 处 理 氢 钟 的 钟 差 数 据 , 将 计算 获得 的 
氨 钟 时 间 尺 度 作 为 Vondrak-Cepek 组 合 滤波 算法 的 
一 组 输入 数据 . 同时 从 实际 光 钟 运行 的 特点 出 发 ， 
模拟 了 多 种 运行 模式 下 的 光 钟 数据 , 将 光 钟 与 参考 
氮 钟 的 模拟 频 差 数据 作为 Vondrak-Cepek 组 合 滤 波 
算法 的 另 一 组 输入 数据 , 综合 计算 获得 光 钟 和 氧 
钟 的 联合 时 间 尺 度 . 最 后 对 不 同 运行 模式 下 的 联 
合 时 间 尺 度 进行 性 能 分 析 , 验证 Vondrak-Cepek 组 
合 滤波 算法 参与 联合 时 间 尺 度 计 算 的 可 行 性 以 及 
评估 联合 时 间 尺 度 相 较 于 协调 世界 时 (Coordinated 
Universal Time, UTC) 的 性 能 情况 . 


2 原子 钟 数 据 模拟 
对 原子 钟 的 数据 模拟 , 我 们 首先 通过 观察 阿 伦 


一 


方差 (Allan Deviation) 相 对 于 观测 区 间 r 的 斜率 , 识 
别 出 原 子 钟 的 主要 噪声 类 型 . 识别 了 主要 噪声 后 佑 
计 其 噪声 参数 值 , 依据 给 定 频 标的 性 能 指标 对 氧 钟 
和 光 钟 进行 模拟 , 利用 实测 数据 和 规格 说 明 书 的 指 
标 进 行 对 比 , 验证 模拟 结果 的 有 效 性 [24. 


2.1 ”主动 型 氧 钟 

中 国 科 学 院 国 家 授时 中 心 (National Time Ser- 
vice Center, NTSC) 是 我 国 唯 一 的 专门 、 全 面 从 事 
时 间 频 率 基 础 研究 和 应 用 研究 的 科研 机 构 , 承担 着 
我 国标 准时 间 UTC (NTSC) 的 产生 、 保 持 和 发 播 
任务 . 利用 我 国标 准时 间 UTC (NTSC) 系 统 中 4 台 
氧 钟 (H1-H4) 的 稳定 度 、 频 率 偏 移 和 频率 漂移 率 进 
行 氨 钟 的 数据 模拟 , 其 中 频率 偏 移 和 频率 漂移 率 参 
考 国 际 权 度 局 (Bureau International des Poids et 
Mesures, BIPM) 网 站 公布 的 数据 . 这 4 台 氢 钟 的 性 
能 指标 如 表 1 所 示 加 ， 


表 1 4 台 和 氢 钟 的 性 能 指标 
Table 1 The performance indicators of four active hydrogen masers 


T Averaging Hydrogen maser Hydrogen maser Hydrogen maser Hydrogen maser 
e 
i time/s H1 H2 H3 HA 
1 1.5 x 1071? 1.6 x 1071? 1.0 x 1071? 9.0 x 10-4 
10 2.5 x 1071 1.8 x 10714 2.5 x 10714 2.5 x 10714 
E 100 6.0 x 10715 4.5 x 10715 7.5 x 10715 6.5 x 10715 
stability uu D oh ii 
i 1000 2.0 x 10 1.5 x 10 2.0 x 10 1.0 x 10 
10000 1.5 x 107 !5 4.5 x 10-16 8.5 x 10-16 9.5 x 10-18 
86400 7.0 x 10716 4.0 x 10-16 7.0 x 10718 6.5 x 10-18 
F 
dod LA z —14 x 107" 5.4 x 10714 4.1 x 1072 74 x 1074 
offset 
Linear 
frequency 一 4 x 10-16 6 x 10718 3.5 x 10718 一 7.0 x 10718 
drift/d 
4E 和 噪声 参数 后 ， 二 不同 的 品 UR EH ite Phase Noise, E 点 状 
确定 了 氢 钟 了 1 的 噪声 参数 后 , 针对 不 同 的 1 白 噪声 (White Phase N WPM), 褐色 点 状 曲 


声 类 型 绘制 出 相应 的 曲线 图 , 如 图 1 所 示 , 蓝 色 星 状 
曲线 代表 毛 钟 H1 的 频率 白 噪声 (White Frequency 
Noise, WFM), 红色 圆圈 曲线 代表 和 氨 钟 H1 的 相位 


线 代 表 氧 钟 H1 的 频率 闪烁 噪声 (Flicker Frequency 
Noise, FFM), 而 黑色 加 号 曲线 代表 氧 钟 HI 的 随机 
游 走 频率 噪声 (Random Walk Frequency Noise, R- 
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WFM). 利用 表 1 中 的 性 能 指标 进行 氨 钟 数据 模拟 ， 模拟 数据 的 合理 性 . 分 别 对 这 4 台 氧 钟 的 主要 噪声 
紫色 曲线 展示 了 和 氧 钟 了 I 的 最 终 模拟 结果 . 这 一 结果 类 型 进行 识别 后 , 估计 的 噪声 参数 值 见 表 2 四. 

与 氨 钟 实际 的 典型 噪声 特性 基本 吻合 , 同时 模拟 计 

算 的 氧 钟 稳定 性 与 实测 结果 基本 一 致 , 从 而 验证 了 
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图 1 氢 钟 1 的 原子 钟 模拟 结果 
Fig.1 Simulation results of hydrogen maser H1 


表 2 4 台 和 氢 钟 的 噪声 参数 
Table 2 The noise parameters of four active hydrogen masers 


Stochastic Hydrogen maser Hydrogen maser Hydrogen maser Hydrogen maser 
component H1 H2 H3 H4 
WPM 1.5 x 107" 1.6 x 107" 1.0 x 107? 9.0 x 107 ^ 
WFM 6.0 x 10 4 4.5 x 107 4 7.5 x 1071 6.5 x 1071 
FFM 6.0 x 10-19 4.0 x 10718 5.0 x 10718 6.5 x 10-15 
RWFM 1.70 x 1071? 9.6 x 10? 1.45 x 107 !5 1.56 x 1075 
Duam T 54 x 1074 4.1 x 10-12 7.4 x 1074 
component 
roni ni 4 x 10719 6 x 10716 3.5 x 10719 —1 x 10-16 
drift/d 
2.2 Xh 图 2 中 模拟 光 钟 与 NIST 光 钟 的 稳定 度 基本 吻 
参考 BIPM 网 站 公布 的 美国 国家 标准 与 计量 研 合 , 同时 在 图 中 可 以 看 出 , 模拟 的 光 钟 噪声 主要 以 
究 院 (National Institute of Standards and Technol- 频率 白 噪 声 为 主 , 也 证 明了 模拟 光 钟 数据 的 有 效 
ogy,，NIST) 的 光 钟 性 能 指标 : 频率 稳定 度 为 3.1 x 性 . 由 于 难以 保持 光 钟 的 整个 复杂 系统 处 于 最 佳 状 
10-16(7/s)-12， 光 钟 几乎 没有 频率 和 频率 漂移 加， 态 , 使 得 光 钟 的 连续 运行 非常 困难 , 因此 考虑 到 光 
其 中 7 代表 了 平均 时 间 , s 表 示 单 位 秒 . 利用 这 台 光 ” 钟 实 际 运行 的 特性 , 我 们 模拟 了 下 列 两 类 情况 . 
钟 的 性 能 指标 进行 模拟 , 模拟 的 光 钟 与 NIST 光 钟 等 间距 模式 : 光 钟 按照 一 个 固定 的 时 间 间 隔 周 


的 稳定 度 对 比 情况 如 图 2 所 示 芭 入 . 


期 性 地 启动 和 停止 . 在 这 种 模式 下 , 固定 周期 的 运 
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行 模式 有 助 于 同步 整个 系统 的 操作 . 

非 等 间距 模式 : 光 钟 启动 之 间 的 时 间 间 隔 可 以 
根据 需要 进行 调整 , 不 遵循 固定 模式 . 这 种 模式 允 
许 更 灵活 的 操作 , 可 以 根据 维护 需求 、 系 统 校 准 或 
其 他 因素 来 调整 停止 时 间 , 更 符合 实际 运行 特点 . 


0 
= — simulated optical clock 


Ú @ NIST optical clock 


Allan Deviation 
全 
4 
/ 
gf 


1 
q 
1 


-18 1 
109 10 


10 


1 4 


10? 10 
Averaging Times/s 


图 2” 模拟 的 光 钟 与 NIST 光 钟 的 频率 稳定 度 对 比 结果 
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Fig.2 Comparison results of frequency stability between 


simulated optical clock and NIST optical clock 


第 一 类 情况 , 等 间距 模式 . 光 钟 运行 期 间 , 模拟 
的 数据 间隔 为 1 h. 
(1) 场 景 1: 理想 情况 下 的 光 钟 数据 , 光 钟 连续 
运行 , 每 小 时 都 有 数据 . 
(2) 场 景 2: 光 钟 每 天 固定 时 段 运 行 1h (运行 
率 为 4.17%); 每 3d 固 定时 段 运行 1h (运行 率 为 
1.3996); 每 7d 固 定时 段 运行 1h (运行 率 为 0.6%). 

(3) 场 景 3: 光 钟 的 运行 率 为 25%, 光 钟 每 天 固 
定时 段 连续 运行 6h; 每 3d 固 定时 段 连续 运行 18h; 
每 7d 固 定时 段 连 续 运 行 42h. 

(4) 场 景 4: 光 钟 的 运行 率 为 50%, 光 钟 每 天 
固定 时 段 连续 运行 12h; 每 3d 固 定时 段 连续 运行 
36h; 每 Yd 固定 时段 连 续 运 行 84h. 

(5) 场 景 5: 光 钟 的 运行 率 为 75%， 光 钟 每 天 
固定 时 段 连 续 运 行 18h; 每 3d 固 定时 段 连 续 运 行 
54h; 每 7d 固 定时 段 连续 运行 126. 

第 二 类 情况 , 非 等 间距 模式 . 光 钟 运行 期 间 , 模 
拟 的 数据 间隔 为 1h. 

(1) 场 景 6: 光 钟 每 天 随机 时 段 运 行 1h (运行 

率 34.1796); 每 3d 随 机 时 段 运行 1h (运行 率 为 


1.3996); 每 7d 随 机 时 段 运行 1h (运行 率 为 0.6%). 

(2) 场 景 7: 光 钟 的 运行 率 为 25%, 光 钟 每 天 随 
机 时 段 连续 运行 6h; 每 3d 随 机 时 段 连 续 运 行 18h; 
每 7dq 随 机 时 段 连续 运行 42h. 

(3) 场 景 8: 光 钟 的 运行 率 为 50%, 光 钟 每 天 
随机 时 段 连续 运行 12h; 每 3d 随 机 时 段 连续 运行 
36h; 每 Yd 随机 时 段 连续 运行 84h. 

(4 场景 9: 光 钟 的 运行 率 为 75%, 光 钟 每 天 
随机 时 段 连续 运行 18h; 每 3d 随 机 时 段 连续 运行 
54h; 每 7d 随 机 时 段 连续 运行 126h. 


3 光 钟 和 和 氢 钟 联合 的 时 间 尺 度 算 法 
氨 钟 具有 优越 的 中 短期 频率 稳定 性 和 连续 稳 
定 运行 的 特点 , 常 被 时 间 基 准 系统 用 为 主 钟 频 率 
i. 氢 钟 的 阿 伦 偏差 在 104 — 105 s 处 达到 最 低 的 频 
率 不 确定 性 , 105 s 至 1 周 左右 , 频率 闪烁 噪声 约 为 2 
一 4 x 10-16 量 级 . 由 于 氢 钟 频 漂 的 存在 , 会 导致 约 
1 周 后 的 频率 稳定 度 变 差 . 
图 3 是 氧 钟 、 钨 喷泉 钟 (Cesium fountain) 和 光 
钟 的 频率 稳定 度 对 比 结果 . 从 图 中 可 以 看 出 , 光 
钟 的 稳定 度 表现 最 优 ， 如果 利用 光 钟 与 氨 钟 的 时 
差 / 频 差 数据 计算 频率 稳定 度 , 最 终 反 映 出 的 是 氧 
钟 的 稳定 度 . 如 果 以 光 钟 作为 参考 , 选用 较 短 时 间 
的 光 钟 与 氧 钟 的 时 差 / 频 差 数据 , 就 可 以 准确 评估 
氧 钟 的 频率 和 频 漂 . 如 果 以 钨 喷泉 钟 为 参考 , 需要 
一 周 以 上 的 时 间 才 能 准确 评估 出 氢 钟 的 频率 和 频 
漂 , 因此 光 钟 参与 时 间 尺 度 计算 具有 很 明显 的 优 
势 


大 多 数 实验 室 的 光 钟 , 每 次 光 钟 停止 运行 到 下 
一 次 开始 运行 的 时 间 间 隔 为 随机 的 , 因此 光 钟 参与 
时 间 尺 度 计 算 带 来 很 明显 的 劣势 . 为 解决 这 一 问题 ， 
部 分 实验 室 采 用 历史 的 光 钟 数据 预测 守 时 型 原子 
钟 的 频率 和 频 漂 , 最 后 利用 守 时 型 原子 钟 计算 出 时 
间 尺 度 . 或 用 内 插 的 方式 修补 光 钟 的 缺失 数据 , 将 
经 过 数据 修补 的 光 钟 作为 一 台 连 续 的 原子 钟 , 与 守 
时 型 原子 钟 联合 计算 获得 时 间 尺 度 . 当 光 钟 运行 率 
较 高 且 缺 失 数据 较 少 时 , 这 两 种 方式 效果 明显 . 但 
光 钟 停止 到 重新 启动 运行 的 时 间 间 隔 较 长 或 缺失 
数据 较 多 时 , 内 插 可 能 会 引入 较 大 的 误差 , 最 终 降 
低 光 钟 参与 计算 的 时 间 尺 度 的 性 能 . 因此 本 文 提 出 
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Vondrak-Cepek 组 合 滤波 算法 解决 光 钟 参与 时 间 尺 
度 计算 的 难点 问题 , 不 需要 考虑 原子 钟 和 光 钟 数据 
的 时 间 间 隔 是 否 均匀 的 , 同时 可 以 消除 原子 钟 的 高 
频 噪声 . 


——Hl -r H2 
m fountain ***9*** Optical clock 


Allan Deviation 
© 


Averaging Time/s 


图 3 Bb. HADER BINDER IR EER ER 


Fig.3 Comparison results of frequency stability among 


hydrogen maser, cesium fountain clock and optical clock 


3.1 Vondrak-Cepek 组 合 滤波 算法 

Vondrak-Cepek 组 合 滤波 算法 相 比 于 Vondrak 
初始 平滑 算法 , 不 但 考虑 了 观测 量 , 还 考虑 了 观测 
量 导数 的 作用 . 其 思路 是 在 曲线 的 光滑 性 与 逼真 度 
之 间 寻 找平 衡 , 该 方法 可 消除 非 均匀 观测 量 中 存在 
的 高 频 噪 声 臣 :1 

假设 有 两 组 观测 数据 , 这 两 组 观测 量 有 不 同 的 
测量 精度 , 其 中 一 组 观测 数据 代表 了 观测 量 , 另 一 
组 观测 数据 代表 了 观测 量 的 一 阶 导 数 , 两 组 观测 是 
不 相关 的 . 定义 了 3 个 变量 : 

(1) 曲 线 的 光滑 性 , ASRR, 如 (1) 式 : 


N-3 
S= 2 (aiyi 
i=l 
其 中 : 


Diyiti + ciyit2 + diYi+3), (1) 


6V (zi42 下 Zi+1) (ZN z 7) 


d (zi — Lit1) (zi — i42) (zi — $443) 
o 6y/(£i42 一 2441) (ZN 一 21) 
(zii = gi) (zi 1 7 Tij 2) (xi 1— Ti 3) : 
PO 6y/(Ti42 一 dd) (ZN 7) 
' (zio = zi) (Zi+2 = Ti41) (i42 一 Tia) ' 
des 6V (isa = Zit1) (ZN ax: xı) (2) 


[tug 一 Li) (Lit3 — Tit1) us — 2) 


yi 为 平滑 值 , x; 为 平滑 值 对 应 的 时 间 序 列 , 入 为 平滑 
值 的 数量 ; 

(2) 平 滑 值 与 观测 量 的 逼真 度 , 用 下 表示 , 如 (3) 
式 所 示 : 


ji ; 2 
F=- p (u) , (3) 
j=1 


中 人 y 为 观测 量 , z) 为 多 对 应 的 时 间 序列 , 每 个 多 有 
不 同 的 权重 p;, n 为 观测 量 的 数量 ; 

(3) 平 滑 值 与 观测 量 一 阶 导数 的 逼真 度 , HER 
示 , in (A) ia: 


认为 观测 量 的 一 阶 导 数 , zu JIA E G h TRI FE AA, 
每 个 冰 有 不 同 的 权重 5, 元 为 观测 量 一 阶 导数 的 数 


B. 采用 最 小 二 乘 方法 最 终 实现 3 个 限制 条 件 最 小 ， 
如 (5) 式 所 示 : 


Q = S +EF +EF = min 

0Q 

Oyi 

这 是 一 组 NN 维 线性 方程 组 , 其 中 e > 0, 5 2 0, 通过 
调整 这 两 个 参数 的 值 , 最 终 实现 @ 值 最 小 . 


3.2 ” 光 钟 参与 计算 的 时 间 尺 度 算法 

利用 单 台 氢 钟 的 性 能 指标 , 模拟 得 到 氧 钟 了 1、 
H2、H3 和 H4 相 对 于 UTC 的 时 差 数 据 , 采用 ALGOS 
算法 计算 获得 氨 钟 时 间 尺 度 TAH 与 氨 钟 HI 的 时 差 
结果 , 该 结果 作为 Vondrak-Cepek 组 合 滤波 算法 的 
一 组 输入 数据 (观测 值 ). 

依据 2.2 节 不 同 场景 下 的 光 钟 模拟 结果 , 可 以 
最 终 计算 获得 光 钟 相对 于 氧 钟 H1 的 时 差 数据 ,并 
将 时 差 数 据 转 化 成 频 差 数 据 , 该 数据 作为 Vondrak- 
Cepek 组 合 滤波 算法 的 另 一 组 输入 数据 (观测 值 的 
一 阶 导数 ). 

该 两 组 数据 是 后 续 联 合 时 间 尺 度 算法 的 数据 
来 源 , 采用 Vondrak-Cepek 组 合 滤 波 算法 , 最 终 计 
算 获 得 氨 钟 和 光 钟 联合 的 时 间 尺 度 TA, 氧 钟 钟 组 
与 光 钟 联合 的 时 间 尺 度 算法 实现 框图 如 图 4 所 示 . 


0,i=1,2 N, (5) 
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A ” 氧 钟 与 光 钟 联合 的 时 间 尺 度 实现 框图 
Fig.4 Block diagram for the time scale implementation combining hydrogen masers and an optical clock 
E I E ps y4 BOE. 
4A 结果 与 分 析 (场景 1, 理想 情况 ) 与 氢 钟 的 频 差 数据 以 及 氧 钟 时 


基于 氧 钟 和 光 钟 的 模拟 比 对 数据 , 利用 上 述 的 
算法 , 针对 不 同 的 光 钟 运行 场景 , 分 别 计算 光 钟 与 
氨 钟 实现 的 组 合 时 间 尺 度 , 最 终 对 联合 时 间 尺 度 进 
4T TERE 4) U2-15]. 

利用 4 台 氢 钟 与 UTC 的 模拟 时 差 数 据 , 采用 
ALGOS 算 法 , 计算 获得 氢 钟 时 间 尺 度 TAn, 如 图 5 
所 示 . 从 图 5 可 以 看 出 , 氢 钟 时 间 尺 度 在 104 一 105 s 
处 达到 噪 底 , 约 为 3 x 10-16. 相 比 于 单 台 原子 钟 , A 
钟 时 间 尺 度 TAH 的 频率 稳定 度 表 现 更 优 . 


102 EN = cde HI 
PIN S S S S S S  LLLOmuS H3 -—wue--H4 
—— TA 
H 
E a 
S -14 、 
glor "Sa. 
3 b A 
a ~ e^ 
~ Do 
s ` P d 
< 101 ~» —'Ü »^? 
Ce 
10-16 i 
10? 10? Br of 108 
Averaging Time/s 
K5 氧 钟 时 间 尺 度 的 频率 稳定 度 比较 


Fig.5 The frequency stability results of time scale based on 


hydrogen masers 


第 一 类 情 ) 
利用 Vondrak-Cepek 组 合 滤 波 算法 , 根据 光 钟 


4.1 


间 尺 度 TAa 与 UTC 的 相差 数据 , 计算 获得 氧 钟 与 光 


钟 联合 的 


间 尺 度 TAa 进 行 分 析 比 较 ， 


时 间 尺 度 TA, 并 将 该 时 间 


结果 如 图 


Scenario 1 
T 7 


尺度 与 氧 钟 时 
6 所 示 . 


UTC-TA/ns 


—— —TA ——TA, 


05 i i , i i i i 
0 20 40 60 80 100 120 140 160 180 
Time/d 
图 6 场景 1 下 计算 的 时 间 尺 度 (TA) 与 UTC 的 时 差 数 据 


Fig.6 Time offsets between time scale TA with UTC for 


Scenario 1 


场景 1 的 光 钟 参与 计算 的 联合 时 间 尺 度 TA 与 


UTC 的 时 间 偏 差 保持 在 


0.1 ns 以 内 , 均 方 根 为 


0.022 ns. 仅 依赖 于 氧 钟 钟 组 计算 的 时 间 尺 度 TAH， 


与 UTC 的 时 间 4 
性 能 明显 差 于 光 钊 


i 差 保 持 在 土 2 ns, 均 方 根 为 0.25 ns, 
参与 计算 的 联合 时 间 尺 度 . 


利用 Vondrak-Cepek 组 合 滤 波 算法 , 根据 光 钟 


(场景 >、 场景 3、 场 景 4 和 场景 5) 与 氨 钟 的 频 差 数据 
以 及 氧 钟 时 间 尺 度 TAH 与 UTC 的 时 差 数 据 , 针对 光 
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MUNI A 与 计算 的 联合 时 
间 尺 度 TA 性 能 , 并 对 这 4 种 情况 进行 分 析 比 较 . 
场景 2 中 ， ,每 天 固定 时 段 运行 1 的 光 钟 参 与 时 
REHA, 产生 的 联合 时 间 尺 度 记 为 TAlinyidi 
每 3d 固 定时 段 运行 1h 的 光 钟 参与 时 间 尺 度 计 
算 , 产生 的 联合 时 间 尺 度 记 为 TAinysa; 每 7d 固 
定时 段 运行 1h 的 光 钟 参与 时 间 尺 度 计算 , 产生 
的 联合 时 间 尺 度 记 为 TAinayvra. 这 3 种 情况 产生 的 
时 间 尺 度 相 比 于 UTC, 其 最 大 绝对 偏差 保持 在 
+0.6 ns 以 内 , 如 图 7 (a) 所 示 . 

场景 3 中 , 每 天 固定 时 段 连续 运行 6h 的 光 钟 
参与 时 间 尺 度 计 算 , 产生 的 联合 时 间 尺 度 记 为 
TAa; 每 3d 固 定时 段 连续 运行 18p 的 光 钟 参 
与 时 间 尺 度 计算 , 产生 的 联合 时 间 尺 度 记 为 


Scenario 2 
; i 


(a) 


UTC-TA/ns 


Ih/1d TAinad 1h/7d 


TA. 
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Time/d 
i Scenario 4 
(c) 
0.5 F 1 
四 
ES] 
t deese edt 
E 0 ] 
I9! 
请 
D 
-0.5 F 
TA mid TA emd TA apa 
E 1 i 1 1 1 i i | 
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Time/d 


TAishy/ad; 每 7d 固 定时 段 连续 运行 42n 的 光 钟 参 
与 时 间 尺 度 计 算 , 产生 的 联合 时 间 尺 度 记 为 
TA4ony7a. 这 3 种 情况 产生 的 时 间 尺 度 相 比 于 UTC， 
其 最 大 绝对 偏差 保持 在 土 0.5 ns 以 内 , 如 图 7 (b) 所 
ZN. 

场景 4 中 , 每 天 固定 时 段 连续 运行 12h 的 光 
钟 参与 时 间 尺 度 计算 ， 产生 的 | 联合 时 间 尺 度 
W TAa; 每 3d 固 定时 段 连续 运行 36h 的 光 
钟 参与 时 间 尺 度 计 算 , 产生 的 联合 时 间 尺 度 
记 为 TAsehyad; 每 7d 固 定时 段 连续 运行 84h 的 光 
钟 参 与 时 间 尺 度 计算 , 产生 的 联合 时 间 尺 度 
记 为 TAganyra. 这 3 种 情况 产生 的 时 间 尺 度 相 比 
于 UTC, 其 最 大 绝对 偏差 保持 在 土 0.54 ns 以 内 , 如 
图 7 (c) 所 示 . 


Scenario 3 
7 r 


(b) 


UTC-TA/ns 


-0.5 - 


TAa; d 18h/3d TA dana 
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图 7 场景 2-5 下 计算 的 时 间 尺 度 与 UTC 的 比较 


Fig.7 Time offsets between time scale TA with UTC for Scenarios 2-5 
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场景 5 中 , 每 天 固定 时 段 连续 运行 18h 的 光 钟 
参与 时 间 尺 度 计算 , 产生 的 联合 时 间 尺 度 记 为 
TAisn/1a; 每 3d 固 定时 段 连续 运行 54n 的 光 钟 参 
与 时 间 尺 度 计 算 , 产生 的 联合 时 间 尺 度 记 为 
TAs54,/34; 每 7d 固 定时 段 连续 运行 126h 的 光 钟 
参与 时 间 尺 度 计 算 , 产生 的 联合 时 间 尺 度 记 为 
TA;»en;za. 这 3 种 情况 产生 的 时 间 尺 度 相 比 于 
UTC, 其 最 大 绝对 偏差 保持 在 土 0.45 ns 以 内 , 如 图 7 
(d) 所 示 . 

同时 , 在 场景 2-5 下 , 光 钟 参与 计算 的 时 间 尺 度 
4232.5 x 10-1 的 较 小 斜率 , 这 主要 是 受 氧 钟 的 漂 
移 影 响 . 


场景 2-5 的 时 间 尺 度 性 能 比较 结果 见 表 3. 从 
表 3 可 以 看 出 , 光 钟 的 运行 率 与 联合 时 间 尺 度 的 性 
能 具有 相关 性 的 , 大 多 数 情 况 下 , 光 钟 的 运行 率 
越 高 , 联合 时 间 尺 度 TA 与 UTC 的 最 大 绝对 偏差 越 
小 , 均 方 根 越 小 . 24369135 fT 838 £7 3 4X 7310096, 
联合 时 间 尺 度 TA 与 UTC 的 最 大 绝对 偏差 在 180 d UJ 
内 , 保持 在 土 0.1 ns, 均 方 根 为 0.022 ns. 而 当 光 钟 的 
运行 率 为 0.6%, 联合 时 间 尺 度 TA 与 UTC 的 最 大 绝 
对 偏差 在 180d 以 内 , 也 能 保持 在 土 0.6 ns, 均 方 根 为 
0.31 ns. 即使 光 钟 的 运行 率 仅 为 0.6%, 它 参 与 计算 
的 时 间 尺 度 性 能 也 明显 优 于 由 纯 氨 钟 参与 计算 的 
时 间 尺 度 TAE. 


表 3 多 种 时 间 尺 度 (Scenario 1-5) 的 性 能 比较 


Table 3 Performance comparison of multiple time scales (Scenarios 1—5) 


Maximum absolute time 


Scenario Optical clock uptime ratio Time scale differences between Root Mean Square 
TA and UTC 
1 100% TAsuM 0.1 ns 0.022 ns 
4.1796 TAin/1a 0.4 ns 0.17 ns 
2 1.39% TAin/3d 0.5ns 0.19ns 
0.696 TA in/7d 0.6ns 0.31 ns 
TAen/ia 0.5ns 0.12 ns 
3 25% TA 1gn/3d 0.4ns 0.10 ns 
TA 42n/7d 0.5ns 0.14ns 
TA12n/1d 0.42 ns 0.10 ns 
4 50% TA 36n/3d 0.45 ns 0.10 ns 
TA gan/7d 0.54 ns 0.12 ns 
TA1gh/1d 0.42 ns 0.10 ns 
5 75% TAsan/3d 0.42 ns 0.10 ns 
TA126n/7d 0.45 ns 0.10 ns 
4.2 ”第 二 类 情况 B. 


利用 Vondrak-Cepek 组 合 滤波 算法 , 根据 光 钟 
(场景 6、 场 景 7、 场 景 8 和 和 场景 9) 与 UTC 的 频 差 数 
据 以 及 氧 钟 时 间 尺 度 TAH 与 UTC 的 时 差 数据 , 针对 
光 钟 不 同 的 4 种 场景 , 分 别 计算 光 钟 参与 计算 的 联 
合 时 间 尺 度 TA 的 性 能 , 并 对 这 4 种 情况 进行 分 析 比 


场景 6 中 , 每 天 随机 时 段 运行 1h 的 光 钟 参与 时 
间 尺 度 计 算 , 产生 的 联合 时 间 尺 度 记 为 TAluia (y: 
每 3d 随 机 时 段 运行 1h 的 光 钟 参与 时 间 尺 度 计算 ， 
产生 的 联合 时 间 尺 度 记 为 TAlnjsa (wy; 每 7d 随 机 时 
段 运行 1h 的 光 钟 参与 时 间 尺 度 计算 , 产生 的 联合 时 
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[8] FRE IRA TA n/a 6). 这 3 种 情况 产生 的 时 间 尺 度 
相 比 于 UTC, 其 最 大 绝对 偏差 保持 在 土 0.7 ns 以 内 ， 
如 图 8 (a) 所 示 . 

场景 7 中 , 每 天 随机 时 段 连续 运行 6h 的 光 钟 
参与 时 间 尺 度 计 算 , 产生 的 联合 时 间 尺 度 记 
为 TA6nyia (wy; 每 3d 随 机 时 段 连续 运行 18h 的 光 
钟 参与 时 间 尺 度 计算 , 产生 的 联合 时 间 尺 度 记 
为 TAigw3q 0); 每 7d 随 机 时 段 连续 运行 42n 的 光 
钟 参 与 时 间 尺 度 计算 , 产生 的 联合 时 间 尺 度 记 
ZJTA428/7a G4). 这 3 种 情况 产生 的 时 间 尺 度 相 比 
TUTO, 其 最 大 绝对 偏差 保持 在 土 0.5 ns 以 内 ， 如 
图 8 (b) 所 示 . 

场景 8 中 , 每 天 随机 时 段 连续 运行 12h 的 光 
钟 参 与 时 间 尺 度 计算 , 产生 的 联合 时 间 尺 度 记 
为 TAiznyia (wy; 每 3d 随 机 时 段 连续 运行 36h 的 光 
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图 8 场景 6-9 下 计算 的 时 间 


钟 参与 时 间 尺 度 计 算 , 产生 的 联合 时 间 尺 度 记 
为 TAs6ny3a (wy; 每 7d 随 机 时 段 连续 运行 84h 的 光 
钟 参与 时 间 尺 度 计 算 , 产生 的 联合 时 间 尺 度 记 
ZJ'TAsanyza (y. 这 3 种 情况 产生 的 时 间 尺 度 相 比 
于 UTC, 其 最 大 绝对 偏差 保持 在 土 0.54 ns 以 内 ,如 
图 8 (c) 所 示 . 

场景 9 中 , 每 天 随机 时 段 连续 运行 18h 的 光 
钟 参与 时 间 尺 度 计 算 , 产生 的 联合 时 间 尺 度 记 
为 TAishjia 0); 每 3d 随 机 时 段 连续 运行 54h 的 光 
钟 参与 时 间 尺 度 计 算 , 产生 的 联合 时 间 尺 度 记 
JJTAsas/3a C); 每 7d 随 机 时 段 连续 运行 126h 的 光 
钟 参与 时 间 尺 度 计 算 , 产生 的 联合 时 间 尺 度 记 
ZJ'TA326n;7a (1). 这 3 种 情况 产生 的 时 间 尺 度 相 比 
TUTO, 其 最 大 绝对 偏差 保持 在 土 0.45 ns 以 内 ,如 
图 8 (d) 所 示 . 
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尺度 TA 与 UTC 的 比较 


Fig.8 Comparison between time scale TA with UTC for Scenarios 6-9 
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同时 , 在 场景 6-9 下 , 光 钟 参与 计算 的 时 间 尺 度 


是 mm 
BELA 


均 有 一 个 较 小 的 斜率 , 3: 
场景 6-9 的 


钟 的 漂移 影响 . 
时 间 尺 度 性 能 比较 结果 见 表 4. 从 
表 4 也 可 以 看 出 , 光 钟 的 运行 率 与 联合 时 间 尺 度 的 


性 能 具有 相关 性 的 , 大 多 数 情 况 下 , 光 钟 的 运行 率 
越 高 , 联合 时 间 尺 度 TA 与 UTC 的 最 大 绝对 偏差 越 


小 , 均 方 根 越 小 . 当 光 钟 在 随机 时 刻 启动 时 , B 


使 光 


钟 的 运行 率 仅 为 0.6%, 联合 时 间 尺 度 TA 与 UTC 的 


最 大 绝对 偏差 在 180 d 以 内 , 也 能 保持 在 3 


上 0.7 ns, 均 


方 根 为 0.2 ns, 明显 优 于 由 


纯 氧 钟 组 成 计 


得 到 的 


时 间 尺 度 TAH. 对 比 表 3 和 表 4, 发 现 表 4 中 的 均 方 


根 优 卫 
受 光 钟 和 氧 钟 数据 的 仿真 误差 以 及 使 
法 中 平滑 参数 的 差异 
Zk EFTIR, Vondrak-Cepek?H A J8 3 $ 


Cepek 组 合 滤波 入 


` Aje 
以 SF 


于 表 3 中 的 均 方 根 , 这 些 差 异 3 


^ 
于 


导致 . 
法 可 


E 要 是 
iVondrak- 


以 


解决 不 连续 的 光 钟 数据 参与 时 间 尺 度 的 问题 , 同时 


能 有 效 提升 光 钟 与 氢 钟 联合 时 间 尺 度 的 1 
间 尺 度 的 性 能 明显 优 于 纯 氯 旬 


表 4 多 种 时 间 尺 度 (Scenario 6-9) 的 性 能 比较 


Table 4 Performance comparison of multiple time scales (Scenarios 6—9) 


FP 钟 组 . 


Maximum absolute time 


生 能 , 其 时 


Scenario Optical clock uptime ratio Time scale differences between Root Mean Square 
TA and UTC 

4.17% TAin/id 0) 0.45 ns 0.10 ns 

6 1.39% TA1n/3d (£) 0.5 ns 0.10 ns 
0.696 TAin/7d (0) 0.7 ns 0.20 ns 

TAonj1a () 0.35 ns 0.10 ns 

7 2596 TA 1gh/3d (r) 0.3 ns 0.10 ns 
"TAa42n/7a (1) 0.5 ns 0.18 ns 

TA12p/1d (r) 0.4 ns 0.10 ns 

8 50% TA 36n/3d (r) 0.4ns 0.10 ns 

TA gap/7d (r) 0.54 ns 0.12 ns 

TA igp/1d (r) 0.42 ns 0.10 ns 

9 75% TA54h/3d (r) 0.42 ns 0.10 ns 
TA126n/7d (r) 0.45 ns 0.10 ns 

5 结论 WEE X6 pb SUPR TUBOS E P AE, 计 


本 文 利 用 Vondrak-Cepek 组 合 滤 波 入 
解决 间歇 运行 的 光 钟 难以 参与 时 间 尺 度 计算 的 难 
题 . 基于 BIPM 发 布 的 实际 毛 钟 和 光 钟 的 性 能 指标 ， 
氧 原子 钟 和 一 台 光 钟 的 数据 . 考虑 到 光 
青 况 , 我 们 模拟 了 光 钟 的 多 


Es 
r1 


模拟 了 4 
钟 实验 室 的 实际 运 


— 


fT 


法 , 尝试 


' 场 景 , 尽 可 能 还 


原 光 钟 的 实际 运行 4 


D. 采用 


HUP 


Vondrak-Cepek 组 合 滤波 算法 , 利 月 


与 UTC 的 
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算 联 合 时 间 尺 度 . 实验 结果 表明 , RI 
光 钟 数据 进行 时 间 尺 度 计算 时 , 在 180d 内 , 理想 情 
况 下 与 UTC 的 偏差 保持 在 4 
场景 下 , 基本 维持 在 3 
时 
Vondrak-Cepek 组 合 滤 


f 


为 基础 


司 尺 度 相 比 UTC, 偏差 


因 


-0.2 ns 以 内 , 而 在 其 人 
EI ns 以 内 . 纯 氧 钟 钟 组 计算 的 


果 持 在 土 1.5 ns 以 内 ， 


波 算法 有 效 提升 了 光 钟 与 氨 


H 


PRATER EREN 


性 和 稳定 性 , 为 未 来 以 光 钟 


建立 更 精确 的 时 间 基 准 黄 定 了 理论 基础 , 


j 模 拟 氧 钟 和 


H 


e 


此 
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Combination of an Optical Clock and Hydrogen Masers for 
Accurate Time Scale Calculation 


ZHAO Shu-hong? DONG Shao-wul23  BAIShan-shan'? QU Lii? ZHANG Ji-hai!? 
ZHANG Yu!? 
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AnBsTRACT The frequency stability and uncertainty of the optical clock have reached a magnitude of 
10775, making it a promising candidate for the next generation of time and frequency standards, and 
it may be used to redefine the international unit "second". The time scale serves as a benchmark for 
accurately and continuously marking the passage of time, being the foundation of high-precision time 
generation. The production of a time scale relies on the continuous and stable operation of atomic clocks; 
however, optical clocks, as laboratory prototype devices, generally cannot operate continuously, thus 
involving optical clocks in time scale calculation presents a challenging issue. The proposed application of 
the Vondrak-Cepek combined filtering algorithm for joint time scale calculations with an optical clock and 
hydrogen masers aims to address the challenges posed by the intermittent operation of an optical clock. 
Initially, the ALGOS algorithm is used on the time difference data of the hydrogen masers to calculate 
and obtain a continuous and stable hydrogen maser clock time scale. Subsequently, the Vondrak-Cepek 
combined filtering algorithm is used to integrate the hydrogen maser time scale with the optical clock data 
to acquire a combined time scale that involves optical clocks in the calculation. Finally, the experimental 
results prove that the Vondrak-Cepek combined filtering algorithm effectively improves the performance 
of the combined time scale of optical clocks and hydrogen masers, with the time deviation of this time 
scale reaching the sub-nanosecond magnitude when compared with Coordinated Universal Time. 


Key words time: time scale, frequency stability; methods: data analysis, Vondrak-Cepek combined 
filtering algorithm; techniques: hydrogen maser, optical clock 
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